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Q Abstract 
O 

The weak convergence theorems of the one- and two-dimensional simple quantum walks, 

o : 

SQW^, d = 1, 2, show a striking contrast to the classical counterparts, the simple random walks, 



00 



SRW^. In the SRW^, distribution of position X(i) of the particle starting from the origin con- 
verges to the Gaussian distribution in the diffusion scaling limit, in which the time scale T and 



spatial scale L both go to infinity with keeping the ratio L/\^T to be finite. On the other hand, in 
the SQW^, the ratio L/T is kept to define the pseudovelocity V(i) = X(i)/t, and then all joint 
moments of the components Vj(t), 1 < j • < d, of V(t) converges in the T = L — > oo limit. The 



limit distributions have novel structures such that they are inverted-bell shaped and the supports 



of them are bounded. In the present paper we claim that these properties of the SQW^- 1 can 
be explained by the theory of relativistic quantum mechanics. We show that the Dirac equation 



with a proper ultraviolet cutoff can provide a quantum walk model in three dimensions, where 
the walker has a four-component qubit. We clarify that the pseudovelocity V(i) of the quantum 
walker, which solves the Dirac equation, is identified with the relativistic velocity. Since the quan- 
tum walker should be a tardyon, not a tachyon, |V(t)| < c, where c is the speed of light, and this 
• i— i . 

restriction (the causality) is the origin of the finiteness of supports of the limit distributions uni- 
versally found in quantum walk models. By reducing the number of components of momentum in 
the Dirac equation, we obtain the limit distributions of pseudovelocities for the lower dimensional 
quantum walks. We show that the obtained limit distributions for the one- and two-dimensional 
systems have the common features with those of SQW^ and SQW^ 2 ). 
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I. INTRODUCTION 



The notion of quantization of random walks has been widely discussed [1-10] and ap- 
plications of quantum walk models have been studied, for example, in information theory 
[11-13] and in solid-state physics [14]. One of the new topics in the study of quantum walks 
is to discuss the relationship between quantum walk models and the relativistic quantum 
mechanics [15-18]. In the present paper, we concentrate on the weak limit theorem of pseu- 
dovelocity of quantum walk, which was first proved by Konno for the one-dimensional simple 
quantum walk [6, 19] and then has been obtained for other models [20-25], and connections 
between the solutions of the Dirac equations and the quantum walk models are reported. 

Simple random walk (SRW) is a discrete-time stochastic process defined on a lattice such 
that particle hopping is allowed only between nearest neighbor sites at each time step. Let 
Z be the set of all integers, Z = {. . . , —2, —1, 0, 1, 2, 3, . . .}. We consider the models on the 
(i-dimensional hypercubic lattices Z d = {x = (x±, . . . ,Xd) '■ Xj G Z, 1 < j < d}. For each 
site in Z d , there are 2d nearest neighbor sites and we put pj be the probability that the 
particle hops to the j-th nearest neighbor site, 1 < j < 2d. The elementary process of the 
SRW^ is determined by a 2rf-component vector p = (p 1: . . . ,p d ) satisfying the conditions 
< Pj < 1, 1 < j < 2d, and Y^f=\Vj = 1- I n other words, if we prepare such a die that it 
has 2d faces and the j-th face appears with probability pj, 1 < j < 2d, in each casting, a 
path of the ci-dimensional SRW (SRW^) can be determined by a sequential random casting 
of this die. 

The ci-dimensional simple quantum walk, SQyV^ d \ is obtained by quantizing the SRW 1 ^. 
In the quantization, the die is replaced by a "quantum die" , which is expressed by a 2d x 2d 
unitary matrix E U(2ci). In the Fourier space, the nearest neighbor hopping is also 
described by using a 2d x 2d unitary matrix representing a shift operator, which is diagonal 



where i = \J— 1 and kj G (ir, ir] denotes the j-th component of the wave number vector k. 

The dynamics of the SQW^ at each time step is then determined by the unitary matrix 



The important consequence of this quantization is that the state of the particle (quantum 
walker) at each time should be represented by a 2<i-component vector-valued wave function, 



S id \k) = diag[e ifcl ,e- 




(1) 




l/ (d) (k) = S {d \k)A (d) , ke( 




(2) 
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to which the 2d x 2d unitary matrix (2) is operated. We assume that there are only one 
quantum walker at the origin O at time t = 0. We let the walker have an initial 2d- 
component qubit q = (q 1: q 2 , • • • , (fed) with qj G C, 1 < j < 2d and J2jti \lj\ 2 — 1- That is, 
in the k-space the initial wave function is independent of k and simply given by 

*(k,0) = * o = t (q 1 ••• Q2 d ), (3) 

where the left-superscript t means the transpose of the matrix. The wave function of the 
walker at time t G N = {0, 1, 2, 3, . . .} is given in the k-space and in the real space R d by 

§(k,t) = [V (d) (k)]'§o (4) 

and 

j pir Ah- _ 

^(x ) t) = n/ 7r ^ k ^(k,t), (5) 

respectively, where k • x = Yfj=\ kjXj. The probability that the quantum walker is observed 
at site x e Z d at time t is given by 

P(x,t) = ||M/(x,t)|| 2 = M/t(x,t)M/(x,t), (6) 

where ^(x, t) = *^(x, t) denotes the Hermitian conjugate of \l/(x, t). In this paper z denotes 
the complex conjugate of z G C. Let X(t) = (Xi(t), X 2 (t), . . . , X d (t)) be the position of the 
quantum walker at time t, whose probability distribution is given by (6). The joint moment 
of the components Xj(t), 1 < j < d, of X(t) is then obtained by [15] 

\j=i / xez d i=i 

d r* dk - ~ d ( 8 \ aj ~ 

= nLf* , * t )ja( i 4) * (M) (7) 

for any aj G N , 1 < j < d. 

The unitarity of the time-evolution operator (2) implies that in principle we are not able 
to find any convergence property of wave function \l/(x, t) nor of the probability distribution 
P(x, t) in the long-time limit t — > oo in the SQW^. It presents a striking contrast to 
the classical stochastic processes, SRW^ , which will generally converge to diffusion particle 
systems in the long-time T — > oo and large-scale L — > oo limit with the diffusion scaling 
L/y/T = const., and probability laws of the diffusion particle systems are described by 



using the Gaussian distribution functions. Konno [6, 19] discovered, however, that in the 
one-dimensional model, SQW^, if we consider the pseudovelocity defined by 



Vi(t) 



(8) 



instead of the position Xi(t) or the usual velocity dXi(t)/dt, it does converge in a weak sense 
such that any moment of V\{t) converges to a moment of a continuous random variable as 
t — > oo, whose distribution is given by a novel probability density function pP-\ That is, 
Konno 's weak convergence theorem is given as follows [6, 19]. Consider the SQW« model 
driven by the quantum die 

a b 



—b a 



eU(2), a,beC, 



\a\ 2 + 



where the initial qubit of the walker is q = (gi, q 2 )- Then for any a\ G N 

, . roc 

lim (V^tr) = / dv 1 v^^ 1 \v 1 ;A^ 1 \q), 



where 



with 



z/ (1) ( Vl ; A™ , q) = ^ { Vl ;\a\)M (1) , q 



^\v i; \a\) 



a\ 2 — vj 



M (1 \A^,q) = 1 



\Qi\ 



l<? 2 | 2 + 



1(|ui| < |a|), 
q^ab + q^ab' 



(9) 



(10) 



(11) 



(12) 



(13) 



Here denotes the indicator function of a condition iv; 1{uj) — 1 if to is satisfied and 
1(uj) = otherwise. Figure 1 shows the density function pP-\vi\ |a|), when a = 1/V2. The 
function /i* 1 ) is now called the Konno density function [10]. 

Recently the weak convergence theorem of the two-component pseudovelocity 

'Xi(t) X 2 (i)\ 



v(0 = (mW)) 



(14) 



, t t 

for the SQW^ 2 \ where X(t) = (X^t), X 2 (t)) denotes the position of the walker, was also 
established [23], where the quantum die is parameterized by p e (0, 1) as 

' ~P 1 - P yjp(l-p) \]p{^ -P) ^ 

1-p 



A^ = 



-P \Jp(1-p) \jp( l ~P) 
ylpQ-p) \jp( l ~P) -(!-p) P 
p(l-p) yjp(l-p) P -(!-p) ) 



e U(4). 



(15) 
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FIG. 1: The one-dimensional density function fjb^(yi; \a\) of limit distribution of pseudovelocity 
for the SQWW (the Konno function), when a = l/y/2. 



For any ai,a 2 G N , 

lim (v^V^t)* 2 ) = dv! r dv 2 v^v^ 2 \v 1: v 2 ;A ( - 2 \q) (16) 



where 



with 



^ 2 \v 1 ,v 2 ;A^,c L ) = ^ 2 \v 1 ,v 2 ;p)M^(v 1 ,v 2 ;A^,c L ) (17) 



V {2) (v u v 2 ;p) 



1X 2 {Vi + V 2 + l)(ui - ^2 + l)(ui + W 2 - l)(^l — U 2 — 1) 

xl(*; 2 /p + *; 2 2 /(l-p)<l). (18) 



The explicit expression of .A/p ■ ) (i>i, t>2; function of v = (v 1: v 2 ), the parameter 

p of the quantum die A^ 2 \ and the four-component initial qubit q = (qi, . . . , g 4 ) is given in 
Sec. III. B in [23]. Figure 2 shows the p — 1/2 case of the density function ^ 2 \ which will 
be a two-dimensional extension of the Konno density function (12). The common feature of 
ji^ and fi^ is that they are inverted-bell shaped on bounded supports, an interval (—a, a) 
and an elliptic region {(i>i,i>2) : vf/p + v 2 /(l — p) < 1}, respectively, as shown in Figs.l and 
2. It is in a big contrast with the Gaussian distributions in d = 1 and d = 2, which are bell 
shaped with unbounded supports R and R 2 describing the diffusion scaling limits of the 
SRW models, the one- and two-dimensional Brownian motions. 

The purpose of the present paper is to clarify the physical meaning of the pseudovelocities 
of quantum walkers, which can have limit distributions in the weak sense even in the quantum 




FIG. 2: The two-dimensional density function fi^ 2 \vi,V2',p) of limit distribution of pseudovelocity 
for the SQW( 2 ), when p = 1/2. 



systems, and the origin of the inverted-bell shape of the limit distribution with bounded 
supports. Our argument is based on the resemblance between the quantum walk models 
and time-evolutions of multi-component wave functions in the relativistic quantum mechanics 
[15-18]. Since ^(k) given by (2) is a unitary matrix, we can assign a Hermitian matrix 
H (d) {k) such that 

yW (k) = e -i^)(k)A (1Q) 

where h is a constant, the Planck constant divided by 2n, and then the time evolution (4) 
of the SQW can be regarded as a solution of the equation 

ih^-$(k,t) = H {d) (k)$(k,t). (20) 

Here time t is now thought to be a continuous variable in R. 

In an earlier paper [15], a relation was reported between the SQW* 1 ) model and the 
Weyl equation, which is given by (20) for the two component wave function, ^(k, t) = 
t (ipi(k, t) ^(k, t)) with the Hamiltonian H^(k) — a • k, where cr = (o^, a 2 , <t 3 ) are the 
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Pauli matrices, 




Note that the Weyl equation and the Dirac equation are the time-evolution equations 
for massless and massive relativistic particles in quantum mechanics, both in the three- 
dimensional continuous space R 3 . So we have to note that the dimensionality d of the 
SQW^ and that of the corresponding relativistic quantum mechanics is different from each 
other in general [15]. Moreover, we have to introduce a proper ultraviolet cutoff in the It- 
space for the quantum mechanics, in order to establish the connection with the SQW^ as 
explained below in detail. 

In Sec. II first we introduce the Dirac equation for a free particle in a usual way [26-28]. 
Then we show that it can provide a model describing motion of a quantum walker having 
a four-component qubit by identifying the time-evolution matrix (2) with the Hamiltonian 
matrix for the Dirac equation explicitly. Then we study the long-time behavior of the 
moments of pseudovelocity for a free Dirac particle starting from the origin. Our initial 
state will be considered as an ideal limit of the highly localized state studied by Bracken, 
Flohr, and Melloy [29]. In this situation we clarify the fact that the pseudovelocity is 
exactly equal to the relativistic velocity. Since the Dirac particle as well as our quantum 
walker should be a tardyon, not a tachyon, its velocity is less than the speed of light c, 



This is the physical origin of the finiteness of supports in the limit distributions of V(t) 
universally found in quantum walk models. In order to ensure the convergence of integrals 
giving joint moments of Vj(t), 1 < j < 3, in t — > oo, we have to introduce an ultraviolet 
cutoff in the theory, which will correspond to the fact that quantum walk models are defined 
on discrete spaces such as, lattices and trees. The introduction of an ultraviolet cutoff here 
will be equivalent with replacement of the delta-function-like initial state by a wave packet 
with a finite size as studied by Strauch [18] for the one-dimensional systems (see Sec. Ill 
C). In Sec. Ill, we modify the Dirac equations by reducing the number of components of 
momentum to describing the quantum walk models in the dimensions lower than d = 3. We 




(22) 
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list up the limit distributions of pseudovelocities of quantum walkers corresponding to these 
free Dirac particles in the dimensions d — 1 and 2 with a proper ultraviolet cutoff. Then we 
compare the results with the limit distributions of the SQW^ with d = 1 and 2 obtained in 
the previous papers [6, 15, 19, 23]. Concluding remarks are given in Sec. IV associated with 
Appendix A. 



II. LIMIT DISTRIBUTION OF DIRAC EQUATION 
A. Dirac equation as a quantum walk model 

The Dirac equation for a free particle with the rest mass m in the three-dimensional real 
space R 3 is given by [26-28] 

i/iJU(x,*) = W*(x,t) (23) 
for the four-component wave function \I/(x, i) with the Hamiltonian operator 

W = 74 c ^4~ + m A > ( 24 ) 
\fc=i k / 

where 7„, z/ = 1,2, 3, 4, are the 4x4 gamma matrices satisfying the algebra 

l^lv + lul^ = ^^ v h, /i, v = 1,2,3,4 (25) 



with the 4x4 unit matrix I 4 . In the present paper, we fix the matrix representations as 
follows, 

/ 



1 -\a k \ 



Ik 



v \<j k J 



, k = 1,2,3, 74 = 



-/ 



(26) 



where a k ,k = 1,2,3, are the Pauli matrices given by (21) and / and are the 2x2 unit 
matrix and the 2x2 zero matrix, respectively. From now on we consider the momentum 
P — Vii Pz) °f the particle, instead of the wave number vector k = (hi, k 2 , k 3 ) of the wave 
function, where the de Broglie relation p = ftk is established. The momentum p provides a 
good quantum number and the solution for a given momentum p is given by a plane wave, 

tf p (x, t) = e ip ^ h - iE ^ h u(p), (27) 

where 



V = |P| = \fpJ+pJ+?3, E{p) = y/ipc)* + (mc 2 )^, (28) 



and u(p) is a four-component vector which satisfies the eigenvalue equation 

E(p)u(p) = W(p)u(p) (29) 

with the Hamiltonian matrix 

H(p) = 74 ^i c "YkPk + mc 2 ^ . (30) 

Then, in the momentum space, the wave function is given as 

^(p,t) = e~ m{p)t/h u(p). (31) 

The Hamiltonian matrix given by (30) can be diagonalized by the Foldy-Wouthuysen-Tani 
transformation [30, 31] as follows: 



H{p) = Uipy'Eip) lA U(p) 



(32) 



with (28) and 
U(P) = 



CTfkPk 



, v E(p) + mc 2 I 4 + i J2 

^2E(p) \ V tiy/E(p)+mc* 

( ■ lE(p) + mc 2 ^ 



JlE(p) 





-CPS 



mc* 



-c(pi— ip2) 



\J E(p)+mc 2 E(p)+mc 2 
-c(pi+ip2) cp 3 



\ ^E(p)+mc 2 y / E(p)+mc 2 



Y E(p)+mc 2 

c(p 1 +ip 2 ) 
\J E(p)+mc 2 

E{p) + mc 2 




c(pi-ip 2 ) \ 
^E(p)+mc 2 



-cp-A 



E(p)+mc 2 





mc* 



(33) 



Then we can see that 



-i«(p)*/» = f i f-i^p); 



n=0 



n! 



ft 



00 i /_i+\ n 
= E-y iT iU{p)-'E{p) 74 C/(p)r 

n=0 \ ^ / 



.n=0 



n! V ft 



f(p) 



f/(p)- 1 diag[e- i£ ^/ n , e -W\ e"^/*, e^^] C/(p) 



(34) 



It implies that the Dirac equation provides a quantum walk model, in which the walker is 
driven by the 4x4 unitary matrix 

V{p) = [/(p^diag 



e -iE( P )/n e -iE(p)/n e iE( P )/n e iE(p)/h 



U(p), 



(35) 



where U(p) is given by (33). 



B. Limit distribution of relativistic velocity 



Let X(£) = (Xi(t),X 2 (t),X 3 (t)) be the position of a free Dirac particle at time t starting 
from the origin at time t — 0. The joint moments of Xj(t), 1 < j < 3 are given by 

'3 \ 3 />oo 



Ux^ty 



n 



2vrft 



(36) 



Let Wj(p) be the j-th row of the matrix U(p) given by (33), 1 < j < 4. We assume the 
initial state as 

§(p,0) = § = t (gi g 2 g 3 ?A (37) 

where qj G C, 1 < j < 4, are independent of p and X^=i |?j| 2 = 1- Then 



$(p,t) = [?7(p)" 1 diag[e" i£(p)/a ,e- i£(p)/a ,e i£(p)/a ,e iS(?,)/?i 
= e-^/^w^^C^p) + w 2 (p)C 2 (p)} 
+e iE ^ h {w 3 (p)C 3 (p) + w 4 (p)C 4 (p)}, 

where C,(p) = [w j (p)]t* , 1 < j < 4. 

For QjGNe {1, 2, • • •}, 1 < j < 3, we see [20] 

3 / d \ aj ~ 



C(p) 



(38) 



-iE( P )t/h 



A % / 



t EU a i +0{t 



Eti° ; 



(39) 



wi(p)Ci(p) + w 2 (p)C 2 (p) 

w 3 (p)C 3 (p) + w 4 (p)C 4 (p) 

Since U(p) is unitary, its row vectors {wj(p)}^ =1 make a set of orthonormal vectors, 

wt(p)w fc (p) = 6 jk , l<j,k<4. (40) 



Then we have 



d 



* f (p.*)n ^ *(M = n 



fdE{p)\ 



dp 



x 



A % / 

EL °i 



I^1(P)| 2 + I^(P)| 2 + (-1)^ 1QJ {I^3(P)| 2 + |C 4 ( P )| 2 } 



(41) 
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The pseudovelocity is defined as 



V(f) = 



, x 1 (t) x 2 (t) x 3 (ty 



t t t j ' 

and we obtain the long-time limit of the joint moments; for aj G N , 1 < j < 3 



iim(nw 



\j=l 



= n 



2nh 



{|Ci(p)| 2 + |C 2 ( P )| 2 } + (-l)^{]C 3 (p)| 2 + |C 4 (p)| 2 } 



4\ / 

Since -E'(p) is the relativistic energy of the particle given by (28), 

- 7 -123 

dPj - c E( P y 3 ~ w - 

We change the variables of integrals in (42) from p^s to v/s by the transformation 

.2 Pj 



(42) 



(43) 



Vj = c 



E(p) 



3 = 1,2,3. 



(44) 



We can show that this map p G R 3 i— > v = (i>i, i>2, 1*3) is one-to-one and the image is an 
interior of a circle with radius c, 



vl + vl + v 2 g < c 2 . 



(45) 



Moreover, we find that the following relations are equivalent with (44), 



mv 



Pj 



=:> J = 1,2,3, 



(46) 



where v — |v|. That is, v is the relativistic velocity of the particle. From this expression 
(46) we can readily calculate the Jacobian associated with the inverse map v i— > p, 

.2 



J = det 



dv-i 



jo 



m 



. d Pk\i<j,k<3 E (P) 5 m 



Associated with the change of variables, (Tj(p)'s are replaced by Cj(v)'s and the integrals 
in (42) are rewritten as 



\j=l 



11 



n 



00 d?;.- 1 



2-kTi J 



3 



{(^(v)! 2 + |c 2 (v)| 2 } + ( -i)£^{|c 3 (v)| 2 + |c 4 (v)| 2 } 



xHv;n(v 2 1+ v 2 2 +vl<c 2 ) 

n r^n<>gL(^L(v;q) 



for a,- G N , 1 < j < 3. Here 



and 



with 



&m = [^) (1 _ v ; /c2)5/2 Hv < c), 



/ mc \ 



v 2 jc 2 fl 2 ' 



-^Dirac(v; q) = 1 + ^ Mj^J^Vj 

i=i 



(47) 



(48) 



(49) 



M^Jsti = 2Re(gig 4 + g 2 g 3 ) A<^(q) = 2Im(gig 4 - g 2 g 3 ), 



(3,2) 



^rKL(q) = 2Re(gig 3 - g2g 4 ), 



(50) 



where Re(z) and lm(z) denote the real part and the imaginary part of z G C, respectively. 



C. Ultraviolet cutoff 



The integrals 



dv-; 



II / II v° j »DL(v)M£L(v; q), K « 2 , a 3 ) e N 3 



(51) 



generally do not converge. In order to obtain finite values of physical quantities the mo- 
mentum space should be restricted, and we set an ultraviolet cutoff by introducing a cutoff 
parameter A > as 

p = |p|<A. (52) 

The range of v = (vi, V2,vs) is then 

A 2 c 2 



v\ + v\ + v\ < 



(mc) 2 + A 2 

We can calculate the normalization constant for finite A as 



3 

n/ 

.7=1 J 



00 dvj ( mc 



c ^2Tih) (l-v 2 /c 2 f/ 2 '\ U " \(mc) 2 + \ 2 J 6tt 2 ^ 3 ' 



rl IV < 



\ 2 C 2 



A 3 



(53) 



(54) 
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Finally the distribution function is determined as 



(3) , .v 3 (mc\ 



v 2 /c 2 fl 2 ' 



V < 



\ 2 C 2 



\ (mc) 2 + A 2 



(55) 




(3) 

FIG. 3: The dependence of ^{ rac on the magnitude v of the pseudovelocity. The solid line for the 



cutoff parameter \/{mc) = 1, and the broken line for A/(mc) = 10, respectively. 



We assume that m is the rest mass of an electron. Figures 3 shows ^i rac (v; A) given by 
(55) as functions of v = \Jv 2 + V2 + v 2 > for the cutoff parameter A/(mc) = 1 (by the 
solid line) and 10 (by the broken line), respectively. They should be compared with Fig.l 
in [29]. If A/(mc) — > oo, that is, the ultraviolet cutoff of the energy E — Ac becomes much 
greater than the rest mass energy mc 2 of the particle, /4)iL c (i>; -M ~* $( v ~~ c )\ the velocity 
of spreading of the Dirac particle outwards from the origin tends to be close to the speed 
of light c [16, 17, 29]. Since quantum walk models are defined on lattices, an ultraviolet 
cutoff is naturally introduced so that A/(mc) = E/{mc 2 ) remains to be small, and thus the 
inverted-bell shaped distributions of pseudovelocities are universally observed. 

The reason that the integrals for the moments (51) do not converge is due to our choice 
of initial state such that a Dirac particle is put at the origin at time t — 0. It will be 
equivalent with setting the initial wave function be a delta function at the origin, which 
is not square-integrable. From this point of view, the introduction of ultraviolet cutoff in 
the present paper can be regarded as a modification of the initial state. In the paper [18], 
Strauch introduced a localization parameter a, which controls the effective size of the initial 
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wave packet describing a Dirac particle and a quantum walker in the models. Dependence 
of the spatial distributions of wave packets in t > on the parameter a and "initial qubit" 
was fully studied for the one-dimensional Dirac equation and the quantum walk models. 
We hope that the present paper will show the importance of his study to understand the 
relationship between the quantum walk models and the relativistic quantum mechanics also 
in the higher spatial-dimensions. 



III. LOWER DIMENSIONAL SYSTEMS 



A. Two-dimensional system 



We consider the Hamiltonian matrix 



w( P ) = [/(p)-'/w + KF 7 4^(p) 



with the two-component momentum p = (pi,P2), where p = |p| = \jp\ + p\ and 

1 / rr, . , n k p k 



U(P) 



2E(p) 



E(p) + mc 2 I A + iJ2 



2E{p) 



k=i yE(p)+ mc 2 
mc 2 

y/E(p) + mc 2 c(pi+ip2) 





c(pi-ip 2 ) \ 
yj E(p)+mc' 2 





-c(pi+ip2) 















) 



(56) 



with E(p) = \J (pc) 2 + (mc 2 ) 2 . Following the similar procedure to that given in Sec. II, we 
will obtain the following result under the ultraviolet cutoff (52) with a parameter A; for 
<Xj e N , j = 1,2, 



\j=l 



A i=l—°° C j=l 



(57) 



with v = J v i + 1>2 and 



^DiLc( v ; A, q) = A*Diac( u ; ^^Dirac^ 0), 



,( 2 ) 



(2) 



(58) 
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where 



(2) / u 1 (mc\ 2 1 



7c 2 ) 



:1 \v < 



A 2 c 2 



\| (mc) 2 + A 2 / ' 



-MgL(v;q) = i + E-^Dd(q) 



(2j) 



with 



-^Dirac(q) = 2Re(gig 4 + g 2 g 3 ), «Mg£.(q) = 2Im( gi g 4 - q 2 q 3 ). 



(59) 
(60) 

(61) 



Figure 4 shows the comparison between (a) of the SQW*- 2 ** given by (17) (see [23] 
for detail), and (b) given by (58). Here we have chosen the parameters as (a) p = 

1/2, q = (1/2, — i/2, —1/2, i/2) and (b) X/(mc) = l,q = (1^/2,1^/2,0,0), respectively. 
In this case, has the symmetry Vi, — v 2 ) = v^(vi, v 2 ), but it is not isotropic 

around the origin. Since we have introduce a simple isotropic cutoff in the momentum 
space (52), i^b-ac is isotropic around the origin. (See a remark in Sec. IV.) Figure 5 shows 
the comparison between (a) with p = 1/2, q = (1/2, i/2, i/2, —1/2) and (b) i^rac w hh 
\/ {mc) = l,q= (-(l + i)/(2 v / 2), -(l + i)/(2v^), (1 + i)/(2 v / 2), (1 -i)/(2 v / 2)), respectively. 
In this case, both distributions show similar anisotropy around the origin. 




J?) 0.5 

Dirac 




FIG. 4: Comparison between (a) of the SQW^ 2 ^ and (b) f^ac- The parameters are chosen as 
(a) p= 1/2, q = (1/2, -i/2, -1/2, i/2) and (b) A/(mc) = l,q= (1/^,1/^,0,0), respectively. 
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FIG. 5: Comparison between (a) of the SQW^ 2 ) and (b) ^bh-ac- The parameters are chosen as 
(a) p = 1/2, q = (1/2, i/2, i/2, -1/2) and (b) A/(mc) = 1, q = (-(l + i)/(2^2), -(l + i)/(2>/2), (1 + 
i)/(2v/2), (1 - i)/(2\/2)), respectively. 



B. One-dimensional system 

By setting the momentum be a scalar pi, we can discuss the one-dimensional system with 
a cutoff parameter A. The obtained limit distribution of the one-component pseudovelocity 
is the following; for cti E N 



lim (Vi(tr) = r -^v^vgUv,; A, q) (62) 



with 



where 



^ffiac^i; A - q) = /^dL^i, a)a<SL(*>i; q), (63) 



/^dLK> A) = J. „ V „ /2 1 I Ui < 



A 2 c 2 



2A (1 - t; 2 /c 2 ) 3 / 2 " " \J (mc) 2 + A 2 1 ' ° ! 1 

-^dLK; q) = 1 + 2^Re(g ig 4 + g 2 g 3 ). (65) 

The result (64) can be compared with Eq.(56) in [17]. 

Figure 6 shows the comparison between (a) of the SQW*- 1 - 1 given by (11) and 
( b ) "Slac S iven b y ( 63 )- Here we have chosen the parameters as (a) a = i 
i^/3/10,q = (l/-s/2,i/v^) and (b) A/(mc) = 3, q = (1/^,1/^,0,0), respectively. 
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Both distributions are symmetric. Figure 7 shows the comparison between (a) with 
a = yj7/10,b = i^/3/10,q = (l/V5,2i/V5) and (b) vgl^ with A/(mc) = 3, q = 
( 1 / vTO, 1 / vTO , 2 / vTO , 2 / vTO) , respectively. Both show similar asymmetry. 




(a) (b) 

FIG. 6: Comparison between (a) i/ 1 ) of the SQW^ and (b) ^DiLc- ^ e parameters are chosen 
as (a) a = iy/7JW,b = i v / 37l0,q = (W2,i/>/2) and (b) A/(mc) = 3,q = (1/^,1/^,0,0), 
respectively. 




FIG. 7: Comparison between (a) i/ 1 ) of the SQW^ and (b) ^rac- The parameters are 
chosen as (a) a = ^JTJ\Q, b = i 73/10, q = (1/V5, 2i/>/5) and (b) A/(mc) = 3,q = 
(1/VlO, 1/VTO, 2/y/W, 2/VlO), respectively. 
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IV. CONCLUDING REMARKS 



In the present paper, we have shown that the novel structures of the limit distributions of 
pseudovelocities in the quantum walk models can be explained by the theory of relativistic 
quantum mechanics. We have studied a very simple system, a single free Dirac particle, 
but we put two special setting when we solve the Dirac equation, in order to clarify the 
connection with the quantum walk models. The first one is the initial condition such that 
the particle is located at the origin at t — 0, and the second one is the introduction of an 
ultraviolet cutoff in the theory. 

Importance of the study of a highly localized state of a free Dirac particle, when we 
consider the connection between the relativistic quantum mechanics and the quantum walk 
models, was clearly demonstrated by Bracken, Flohr, and Mellow [29], by Strauch [16], and 
by Bracken, Ellinas, and Smyrnakis [17]. By observing the time evolution of the wave- 
packet-type solutions of the Dirac equation, they found that such a highly localized initial 
state of a free Dirac particle leads to a rapid expansion of the distribution of position of 
the particle outwards from the origin, at speeds close to c. Our initial condition will be 
regarded as an ideal limit of their setting, in which the width of the wave-packet goes to 
zero. If we study the long-time limit of the pseudovelocity V(t) of the free Dirac particle, 
this phenomenon gives a singular distribution such that only on the surface of the sphere 
with the radius c we have the delta measure, and nothing on any other points in the V-space. 
Since the quantum walk models studied so far are defined on discrete spaces, such as lattices 
and trees, we have to introduce an ultraviolet cutoff in the quantum mechanics. Then the 
limit distributions of V(t) are moderated and we have obtained the inverted-bell shaped 
distribution functions. Introduction of an ultraviolet cutoff A < oo in the momentum space 
(52) will be equivalent with the introduction of an effective size a oc h/X > for the spatial 
initial state of a quantum particle/walker. For the one-dimensional models, systematic 
study of the dependence of solutions on the parameter a was reported by Strauch [18] for the 
discrete- and the continuous-time quantum walks and for the relativistic and non-relativistic 
quantum mechanics. 

In the present paper, we have considered an isotropic cutoff in the momentum space by 
introducing a single parameter X; p = |p| < A. We expect that the variety of the quantum 
walk models depending on lattice structure on which a quantum walker exists and choice of 
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a quantum die represented by a unitary matrix will be partially realized by changing how to 
introduce cutoff in the theory. The present isotropic cutoff is too simple and thus we have 
some differences in symmetry of distribution functions between the graphs (a) of for the 
SQW^ and (b) of 

^dILc obtained from the Dirac equations in Figs. 4 and 5 in the d — 2 

case. 

In the previous paper on the SQW^ 2 ^ [23] an interesting phenomenon called the localization 
of quantum walker at the starting point was reported (see also [21, 22]). In the present 
models derived from the Dirac equation, however, the localization phenomenon will not be 
expected. 

In the present study, we have derived the limit distributions for the d — 1 and d = 

2 systems from the result for the original Dirac equation in the three dimensions. The 

dependence of the results on the dimensionality d will be summarized as the following, 

(d) , \ \ f mc\ d 1 (v [ ( mc\ 2 } . . 

^L(v,X) = const. ^— j (1 ^ 2/c2)(rf+2)/2 l y~ < [1+ {— ) ) J • (66) 

Note that, if we define the fifth gamma matrix by 75 = 71727374 and consider the Hamilto- 
nian operator of the form 

n = 74 (s chlk ^ + chl5 ^ + m ° 2 ) ' (67) 

we can discuss the quantum walk in the four-dimensional space with the four-component 
qubit. Also in this case, we can obtain the limit distribution, which has the form (66) with 
d = 4 as shown in Appendix A. 
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APPENDIX A: FOUR-DIMENSIONAL SYSTEM 

The diagonalization of the Hamiltonian matrix corresponding to (67) can be done in the 
four-component momentum space p = (pi,P2,P3,Ps) as follows, 

W(p) = U( P y l E(p) l4 U(p) (Al) 
19 



with 

C^fkPk 



U(p) = -jJ= l^E( P )+mc*I 4 + i £ ^ , (A2) 



where £(p) = ^(p? + p| + p§ + p!)c 2 + (mc 2 ) 2 . Let V(f) = (Vi(t), . . . , V 4 (t)) be the four- 
component pseudovelocity. When we introduce an isotropic cutoff in the model (52), the 
weak limit theorem is obtained as follows; for any aj G N , 1 < j < 4, 



with 



where 



with 



( n v s (t) ai ) = n / n ^Uv; a, q ) (a 3 ) 



\j=l / A j=l—°" ^ j=l 



^blc(v; A, q) = ^L(v, A)-MgL(v; q) (A4) 



a*dL(«; A) = 5r ( ^ I r, — itt^ 1 1 u < 



A 2 c 2 



(A5) 



7T 2 V A J (f - w 2 /c 2 ) 3 ^ \| (mc) 2 + A 2 

-Md1c(v; q) = 1 + E M&Jdvi (A6) 
i=i 



•^DimUq) = 2Re(gig 4 + g 2 g 3 ), ^^(q) = 2Im(g!g 4 - q 2 q 3 ), 

Mggl(q) = 2Re( qi q 3 - q 2 q 4 ), A<8&(q) = 2Im(g ig3 + (fcfc). (A7) 
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